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Abstract. There has been great progress in recent years in discovering star forming galaxies at 
high redshifts (z > 5), close to the epoch of reionization of the intergalactic medium (IGM). The 
WFC3 and ACS cameras on the Hubble Space Telescope have enabled Lyman break galaxies to 
be robustly identified, but the UV luminosity function and star formation rate density of this 
population at z = 6 — 8 seems to be much lower than at z = 2 — 4. High escape fractions and a 
large contribution from faint galaxies below our current detection limits would be required for 
star- forming galaxies to reionize the Universe. We have also found that these galaxies have blue 
rest-frame UV colours, which might indicate lower dust extinction at z > 5. There has been 
some spectroscopic confirmation of these Lyman break galaxies through Lyman-a emission, but 
the fraction of galaxies where we see this line drops at z > 7, perhaps due to the onset of the 
Gunn-Peterson effect (where the IGM is opaque to Lyman-a). 
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1. Introduction 

In the past decade, the quest to observe the most distant galaxies in the Universe has 
rapidly expanded to the point where the discovery of z ~ 6 star-forming galaxies has 
now become routine. Deep imaging surveys with the Hubble Space Telescope (HST) and 
large ground based telescopes have revealed hundreds of galaxies at z ~ 6 (Bunker et 
al. 2004, Bouwens et al. 2006, 2007). These searches typically rely on the Lyman break 
galaxy (LBG) technique pioneered by Steidel and collaborators to identify star-forming 
galaxies at z « 3 — 4 (Steidel et al. 1996, 1999). Narrow-band searches have also proved 
successful in isolating the Lyman-a emission line in redshift slices between z = 3 and 
z = 7 (e.g., Ouchi et al. 2008; Ota et al. 2008). Rapid follow-up of Gamma Ray Bursts 
has also detected objects at z ~ 6, and most recently one spectroscopically confirmed at 
z = 8.2 (Tanvir et al. 2009). 

Parallel to these developments in identifying high redshift objects has been the discov- 
ery of the onset of the Gunn-Peterson (1965) effect. This is the near-total absorption of 
the continuum flux shortward of Lyman-a in sources at z > 6.3 due to the intergalac- 
tic medium (IGM) having a much larger neutral fraction at high redshift. The Gunn- 
Peterson trough was first discovered in the spectra of SDSS quasars (Becker et al. 2001, 
Fan et al. 2001, 2006). This defines the end of the reionization epoch, when the Universe 
transitioned from a neutral IGM. Latest results from WMAP indicate the mid-point of 
reionization may have occurred at z ~ 11 (Dunkley et al. 2009). The source of necessary 
ionizing photons remains an open question: the number density of high redshift quasars 
is insufficient at z > 6 to achieve this (Fan et al. 2001, Dijkstra et al. 2004). Star- forming 
galaxies at high redshift are another potential driver of reionization, but we must first 
determine their rest-frame UV luminosity density to assess whether they are plausible 
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sources; the escape fraction of ionizing photons from these galaxies, along with the slope 
of their UV spectra, are other important and poorly-constrained factors in determining 
whether star formation is responsible for the ionization of the IGM at high redshift. 

2. The Star Formation Rate at High Redshift 

Early results on the star formation rate density at z « 6 were conflicting, with some 
groups claiming little to no evolution to z ~ 3 (Giavalisco et al. 2004, Bouwens et al. 2003) 
while other work suggested that the star formation rate density at z ~ 6 was significantly 
lower than in the well-studied LBGs at z = 3 — 4 (Stanway et al. 2003). The consensus 
which has now emerged from later studies is that the abundance of luminous galaxies is 
substantially less at z « 6 than at z « 3 (Stanway et al. 2003, Bunker et al. 2004, Bouwens 
et al. 2006, Yoshida et al. 2006, McLure et al. 2009). If this trend continues to fainter 
systems and higher redshifts, then it may prove challenging for star- forming galaxies to 
provide the UV flux needed to fulfill reionisation of the intergalactic medium at (e.g. 
Bunker et al. 2004). Importantly, analysis of the faint-end of the luminosity function 
at high redshift has revealed that feeble galaxies contribute an increased fraction of the 
total UV luminosity, suggesting that the bulk of star formation (and hence reionizing 
photons) likely come from lower luminosity galaxies not yet adequately probed even in 
deep surveys. 

Until recently, extending this work to the z ~ 7 universe has been stunted by small 
survey areas (from space) and low sensitivity (from the ground). There is evidence of 
old stellar populations in z ~ 4 — 6 galaxies from measurements of the Balmer break in 
Spitzer/IRAC imaging (Eyles et al. 2005,2007; Stark et al. 2007, 2009), which must have 
formed at higher redshift (see Fig. 3). However, the age-dating and mass-determination of 
these stellar populations has many uncertainties, so searching directly for star formation 
at redshifts z ^ 7 is critical to measure the evolution of the star formation rate density, 
and address the role of galaxies in reionizing the universe. 

The large field of view and enhanced sensitivity of the recently-installed Wide Field 
Camera 3 (WFC3) on HST has made great progress in identifying larger samples of 
z > 6 Lyman break galaxies, as it covers an area 6.5 times that of the previous-generation 
NICMOS NIC3 camera, and has better spatial sampling, better sensitivity and a filter set 
better tuned to identifying high-redshift candidates through their colours. In Bunker et 
al. (2010) we presented an analysis of the recently-obtained WFC3 near-infrared images 
of the Hubble Ultra Deep Field (UDF). We have previously used the i'-band and z' - 
band ACS images to identify LBGs at z ~ 6 through their large %' — z' colours (the 
i'-drops, Bunker et al. 2004). Using this ACS z'-band image in conjunction with the new 
WFC3 F-band we could search for galaxies at z ~ 7, with a spectral break between 
these two filters (the z'-drops). We also anaylsed the J and i^-band WFC3 images to 
eliminate low-redshift contaminants of the z'-drop selection through their near-infrared 
colours (Fig. 1), and also to look for Lyman-break galaxies at higher redshifts (the Y- 
and J-drops at z ~ 8 & z ~ 10). 

3. The Role of Star Forming Galaxies in Reionization 

We can use the observed F-band magnitudes of objects in the z'-drop sample to 
estimate their star formation rate from the rest-frame UV luminosity density. In the ab- 
sence of dust obscuration, the relation between the flux density in the rest-UV around 
& 1500 A and the star formation rate (SFR in M yr _1 ) is given by Ljjv = 8 x 
lO^SFRergss' 1 Hz' 1 from Madau, Pozzetti & Dickinson (1998) for a Salpeter (1955) 
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stellar initial mass function (IMF) with 0.1 M < M* < 125 M . This is comparable 
to the relation derived from the models of Leitherer & Heckman (1995) and Kennicutt 
(1998). However, if a Scalo (1986) IMF is used, the inferred star formation rates will 
be a factor of ~ 2.5 higher for a similar mass range. The redshift range is not surveyed 
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Figure 1. Left: The z'-band dropouts at z ~ 7 identified in the WFC3 and ACS imaging of 
the Hubble Ultra Deep Field (Bunker et al. 2010). The top candidate, zDO, is actually an inter- 
mediate-redshift supernova (rather than a high redshift Lyman break galaxy) which occurred 
during the more recent WFC 3 imaging from 2009 but was absent in the 2003 ACS images. 
Top Right: The WFC 3 and ACS filter set on HST used, along with the spectral template for 
a z ~ 7 star-forming galaxy, and two potential interloper populations: a low mass star in our 
own Galaxy, and an evolved galaxy at z ~ 2.5 (from Wilkins et al. 2011a). 

Bottom Right: The use of a two-colour diagram to reject the low-redshift interlopers and select 
candidate z ~ 8 y-drop galaxies (Lorenzoni et al. 2011). The dotted line is the locus of evolved 
galaxy colours as a function of redshift; the lower points are stellar colours at zero redshift, and 
the near- vertical lines denote star- forming galaxies with spectral slopes /3 = — 3 (left) and f3 — 
(right). The points with error bars are candidate high redshift galaxies in the HUDF09 field 
from Lorenzoni et al. (2011). 
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with uniform sensitivity to UV luminosity (and hence star formation rate). We have a 
strong luminosity bias towards the lower end of the redshift range, due to the effects of 
increasing luminosity distance with redshift and also the Lyman-a break obscuring an 
increasing fraction of the filter bandwidth. Using an "effective volume" approach (e.g., 
Steidel et al. 1999) accounts for this luminosity bias, we infer a total star formation rate 
density at z « 7 of 0.0034 M yr" 1 Mpc -3 , integrating down to 0.2 L^ v of the L* value 
at z = 3 (i.e. Mjjy = —19.35) and assuming a UV spectral slope of f3 = —2.0. If in- 
stead we integrate down to 0.1L* Z=3 (0.2L* =6 ) and the total star formation rate density 
is 0.004 M yr _1 Mpc~ 3 . These star formation rate densities are a factor of ~ 10 lower 
than at z ~ 3 — 4, and even a factor of 1.5 — 3 below that at z ~ 6 (Bunker et al. 2004; 
Bouwens et al. 2006). 

We can compare our measured UV luminosity density at z ~ 7 (quoted above as 
a corresponding star formation rate) with that required to ionize the Universe at this 
redshift. Madau, Haardt & Rees (1999) give the density of star formation required for 
reionization (assuming the same Salpeter IMF as used in this paper): 

0.005 Mpyr^Mpc- 3 f l + z \ 3 f Q b h 2 70 \ 2 (C 



PSFR ~ u c \rr) Komm) U 

We have updated equation 27 of Madau, Haardt & Rees (1999) for a more recent con- 
cordance cosmology estimate of the baryon density from Spergel et al. (2003). C is the 
clumping factor of neutral hydrogen, C = (phj) (Ph/)~ 2 - Early simulations suggested 
C « 30 (Gnedin & Ostriker 1997), but more recent work including the effects of re- 
heating implies a lower concentration factor of C ~ 5 (Pawlik et al. 2009). The escape 
fraction of ionizing photons (f eS c) for high-redshift galaxies is highly uncertain (e.g., Stei- 
del, Pettini & Adelberger 2001, Shapley et al. 2006), but even if we take f esc — 1 (no 
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Figure 2. Left: the rest-frame UV luminosity function of star-forming galaxies at z ~ 7 is 
dramatically lower than at z ~ 2, and even than z ~ 6 (from Wilkins et al. 2011a). 
Top Right: The evolution in the rest-UV colours of bright Lyman break galaxies; at higher 
redshifts, the galaxies are much more blue (from Wilkins et al. 2011b). 

Bottom Right: Various factors can affect the rest-UV spectral slopes of star-forming galaxies, 
such as IMF, metallicity, star formation history and dust (from Wilkins et al. 2011b). The blue 
colours at high redshift might imply these early galaxies were much less dust reddened than the 
corresponding lower-redshift populations. 
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absorption by H i) and a very low clumping factor, this estimate of the star formation 
density required is a factor of 1.5-2 higher than our measured star formation density at 
z ~ 7 from z'-drop galaxies in the UDF. For faint end slopes of a — 1.8 — 1.6 galaxies 
with L > 0.2 L* account for 24 — 44% of the total luminosity (if there is no low- luminosity 
cut-off for the Schechter function), so even with a steep faint-end slope at z ~ 7 we still 
fall short of the required density of Lyman continuum photons required to reionize the 
Universe, unless the escape fraction is implausibly high (f esc > 0.5) and/or the faint end 
slope is a < —1.9 (much steeper than observed at z = — 6). At z ~ 8, the situation 
is even more extreme (Lorenzoni et al. 2011), as the luminosity function is even lower 
(Fig. 4). 

However, the assumption of a solar metallicity Salpeter IMF may be flawed: the colours 
of z ~ 6 z'-band drop-outs are very blue (Stanway, McMahon & Bunker 2005), with 
/3 < — 2, and the new WFC3 J- and i^-band images show that the z ~ 7 z'-drops also 
have blue colours on average (Figure 2). A slope of j3 < — 2 is bluer that for continuous 
star formation with a Salpeter IMF, even if there is no dust reddening. Such blue slopes 
could be explained through low metallicity, or a top-heavy IMF, which can produce 
between 3 and 10 times as many ionizing photons for the same 1500 A UV luminosity 
(Schaerer 2003 - see also Stiavelli, Fall & Panagia 2004). Alternatively, we may be seeing 
galaxies at the onset of star formation, or with a rising star formation rate (Verma et al. 
2007), which would also lead us to underestimate the true star formation rate from the 
rest-UV luminosity. We explore the implications of the blue UV spectral slopes in z ^ 6 
galaxies in Wilkins et al. (2011b). 

We have also targeted many of our high redshift Lyman break galaxies with spec- 
troscopy from the ground, using Keck/DEIMOS (Bunker et al. 2003), Gemini/GMOS 
(Stanway et al. 2004) and most recently VLT/XSHOOTER (Bunker et al. 2012, Fig. 5). 
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Figure 3. Left: the spectral energy distribution from Spitzer/IRAC and HST data for a z ~ 6 
star- forming galaxy (Eyles et al. 2005). If we fit two stellar populations (ongoing star formation 
and an older burst), we find that there are stars as old as 450Myr at z ~ 6 (within 900Myr of the 
Big Bang). The pushes the formation epoch back to z'simlO, around the time to reionization. 
Right: The assembly of stellar mass at high redshift, derived from Spitzer measurements of 
Lyman break galaxies (Eyles et al. 2007). Our high redshift points are lower limits because 
the Lyman break selection only picks out actively star-forming galaxies (and not post-burst 
systems) . 
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While we have seen Lyman-a emission is several z « 6 galaxies, the frequency seems to 
drop at z > 7, perhaps due to absorption of this line in the increasingly-neutral IGM. 



4. Conclusions 

The availability of WFC3 and ACS on the Hubble Space Telescope has made the 
discovery of Lyman break galaxies at z > 5 possible. From these galaxies, we have 
shown that the star formation rate density declines at z > 5. Based on their rest-frame 
UV luminosity function, it appears that a high escape fraction for Lyman continuum 
photons and a huge contribution from faint sources well below current detection limits 
would be needed for star forming galaxies to keep the Universe ionised at these redshifts. 
We have been able to measure stellar masses and ages from Spitzer in some instances. We 
also have discovered that the high redshift Lyman break galaxies have much bluer rest- 
frame UV colours than the z ~ 4 population, perhaps indicative of less dust obscuration 
at higher redshifts. Spectroscopic follow-up is ongoing, but we have confirmed Lyman-a 
emission from several z ~ 6 galaxies, although the success rate drops at z > 7, perhaps 
due to the onset of Gunn-Peterson absorption. 
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Figure 4. Left: the evolution of the star formation rate density (the "Madau-Lilly" diagram) 
based on Lyman break galaxies. Our work at z ~ 6 (Bunker et al. 2004) indicates a much lower 
level of activity than at z ~ 2 — 4. 

Right: the ionising flux produced (y-axis) by the observed luminosity function of star-forming 
galaxies at z ~ 8, integrated down to a faint limiting magnitude (plotted on the x-axis). Where 
the integrated Schechter function curves cross the horizontal lines, there is sufficient ionising 
photons to keep the Universe ionised. As can be seen, a high escape fraction is required, and 
most of the photons would have to come from faint galaxies well below our observational limit 
(the shaded grey box in the lower left). Figure from Lorenzoni et al. (2011). 
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z > 7, the spectroscopy has been less successful, perhaps because the resonant Lyman-a line 
is absorbed by optically-thick neutral hydrogen in the Inter- Galactic Medium during the Gun- 
n-Peterson era. 

Right: we have been unable to confirm the claimed spectroscopic detection of Lyman-a at 
z = 8.55 by Lehnert et al. (2010) with SINFONI in the HUDF F-band drop-out YD3 from 
Bunker et al. (2010). Our 5-hour XSHOOTER spectrum is shown at the top, and Gaussian-s- 
moothed below. The expected position of the line emission is denoted by a circle. The bottom 
two panels have fake sources of the same flux as the Lehnert et al. emission line added in, with 
two different velocity widths. We should have detected the line at « 4 a level if it was real (from 
Bunker et al. 2012). 
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